Climate instability on tidally locked exoplanets 
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ABSTRACT 

Feedbacks that can destabihze the chmates of synchronously-rotating rocky planets may 
arise on planets with strong day-night surface temperature contrasts. Earth- like hab- 
itable planets maintain stable surface liquid water over geological time. This requires 
equilibrium between the temperature-dependent rate of greenhouse-gas consumption 
by weathering, and greenhouse- gas rcsupply by other processes. Detected small-radius 
exoplanets, and anticipated M-dwarf habitable-zone rocky planets, are expected to be 
in synchronous rotation (tidally locked). In this paper we investigate two hypothetical 
feedbacks that can destabilize climate on planets in synchronous rotation. (1) If small 
changes in pressure alter the temperature distribution across a planet's surface such 
that the weathering rate goes up when the pressure goes down, a runaway positive feed- 
back occurs involving increasing weathering rate near the substellar point, decreasing 
pressure, and increasing substellar surface temperature. We call this feedback enhanced 
substellar weathering instability (ESWI). (2) When decreases in pressure increase the 
fraction of surface area above the melting point (through reduced advective cooling of 
the substellar point), and the corresponding increase in volume of liquid causes net 
dissolution of the atmosphere, a further decrease in pressure will occur. This substellar 
dissolution feedback (SDF) can also cause a runaway climate shift. We use an idealized 
energy balance model to map out the conditions under which these instabilities may 
occur. In this simplified model, the weathering runaway can shrink the habitable zone, 
and cause geologically rapid 10^-fold atmospheric pressure shifts within the habitable 
zone. Mars may have undergone a weathering runaway in the past. Substellar dissolu- 
tion is usually a negative feedback or weak positive feedback on changes in atmospheric 
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pressure. It can only cause runaway changes for small, deep oceans and highly soluble 
atmospheric gases. Both instabilities are suppressed if the atmosphere has a high ra- 
diative efficiency. Our results are most relevant for atmospheres that are thin, have low 
greenhouse-gas radiative efficiency, and have a principal greenhouse gas that is also the 
main constituent of the atmosphere. These results identify a new pathway by which 
habitable-zone planets can undergo rapid climate shifts and become uninhabitable. 

13 Subject headings: planetary systems - planets and satellites: general- stars: individual (Kepler- 

14 10, CoRoT-7, GJ1214, 55 Cnc, Kepler-9, Kepler-11) 



1. Introduction 



16 ipxoplanet research is driven in part by the hope of finding habitable planets beyond Earth (Ex- 



oplanet Community Report 2009). Demonstrably habitable exoplanets maintain surface liquid 
water over geological time. Earth's long-term climate stability is believed to be maintained by 
a negative feedback between control of surface temperature by partial pressure of CO2 (PCO2), 



20 and temperature-dependent mineral weathering reactions that reduce pC02 (Walker et al. 1981). 



21 There is increasing evidence that this mechanism does, in fact, operate on Earth (Cohen et al. 



(Walker et al.|1981). 


Earth ( 


Cohen et al. 



22 (2004) Zeebe k Caldeira (2008), but see also Edmond k Huh (2003)). The Circumstellar Hab- 



itable Zone hypothesis (Kasting et al. 1993) extends this stabilizing feedback to rocky planets in 
general, between top-of-atmosphere flux limits set by the runaway greenhouse (upper limit) and 
condensation of thick CO2 atmospheres (lower limit). H2-H2 collision-induced opacity can extend 
the habitable zone further out, in theory ( [Pierrehumbert k Gaidos 2011, Wordsworth||20lT ) . Cur- 



rently the best prospects for finding stable surface liquid water orbit M stars (Tarter et al. 2007 
Deming k Seager||2009 ) . Planets in M-dwarf habitable zones are close enough to their star for tidal 
despinning and synchronous rotation (Murray k Dermott (1999), Chapter 5). Nearby M-dwarfs 
are the targets of several ongoing and proposed planet searches. Rocky exoplanets in hot orbits 



31 have recently been confirmed (Leger et al. 2011, Batalha et al. 2011 Winn et al. 2011), and are 



32 presumably in synchronous rotation. But does the habitable zone concept hold water for tidally 

33 locked planets? 

34 In this paper, we highlight two closely-related feedbacks which could cause climate destabiliza- 

35 tion on planets with and low-opacity atmospheres and atmospheres that do not have large optical 

36 depth. Both feedbacks require surface temperatures near the substellar point to be significantly 

37 higher than the planet-average surface temperature. 



38 
39 
40 
41 



The enhanced substellar weathering instability (ESWI) flows out of the same strong temper- 
ature dependence of silicate weathering that makes it possible for carbonate-silicate feedback 



to stabilize Earth's climate (Walker et al. 1981). Weathering and hence CO2 drawdown rate 



increases rapidly with increasing temperature. Weathering also increases with rainfall, which 
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increases with temperature (O'Gorman & Schneider 2008 Pierrehumbert 2002 Schneider et 



ar]|2010 ). Therefore, the global CO2 loss rate depends heavily on the maximum temperature. 
On a synchronously-rotating planet where ATg is high, most of the weathering occurs near 
the substellar point. Suppose weathering is initially adjusted to match net supply of green- 
house gases by other processes (e.g., volcanic degassing). Then consider a small increase in 
atmospheric pressure. Average temperature must increase, unless there is an antigreenhouse 
effect. Normally, this would lead to an increase in weathering. However, on a planet with 
a strong day- night temperature contrast, most of the weathering occurs near the substellar 
point. An increase in atmospheric pressure can decrease temperature around the substellar 
point, if increased advection of heat away from the hotspot by winds outweighs any increase 
in greenhouse forcing. Because this substellar area is cooling, and most of the weathering is 
around the substellar point, the planet- averaged weathering rate declines. Volcanic supply 
of greenhouse gases now outpaces removal by weathering, and a further increase in pressure 
occurs. This instability can lead to very strong greenhouse forcing, and may trigger a moist 



runaway greenhouse (Kasting 1988). Conversely, a small decrease in pressure from the un- 



stable equilibrium can lead to atmospheric collapse. ESWI requires that weathering is an 
important sink for the major climate-controlling greenhouse gas, which is also the dominant 
atmospheric constituent. It also requires that the atmosphere is important in setting the 
mean surface temperature. 

Substellar dissolution feedback (SDF) supposes an increasing gradient in surface temperature 
on an initially frozen planet, which allows a liquid phase to form (or be uncovered) around 
the substellar point. Some atmosphere dissolves in the new liquid phase. Positive feedback 
is possible if the decrease in atmospheric pressure due to dissolution raises the temperature 
around the substellar point, increasing the fraction of the planet's surface area above the 
melting point. (We assume for the moment that P exceeds the triple point). In order for the 
mass of atmosphere sequestered in the pond to increase with decreasing pressure, increasing 
pond volume must outcompete both Henry's-law decrease in gas dissolved per unit volume, 
and the decrease in gas solubility with increasing temperature. For example, suppose c oc P, 
where c is the concentration of gas in the pond and P is the atmospheric pressure, and 
V oc with V the ocean volume. Then n > 1 is sufficient for positive feedback, and n > 2 
is sufficient for runaway. So long as the runaway condition is satisfied, the area of liquid 
stability will continue to expand: a pond becomes an ocean, drawing down the atmosphere. 
As with the ESWI, the key is that substellar temperature increases as pressure decreases. 
Runaway SDF implies a climate bistability for a given inventory of volatile substance. One 
equilibrium has all of the volatile in the atmosphere. The other equilibrium has most of the 
volatile substance sequestered in a regional ocean and a little in the atmosphere, with the 
ocean prevented from completely freezing over by the steep temperature gradient that the 



thin atmosphere enables. A similar hysteresis was proposed for ancient Titan by Lorenz et 



al. (1997). Runaway SDF is separate from the feedback between retreating ice cover and 



increased absorption of sunlight (ice-albedo feedback; Roe &: Baker (2010)), although both 
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82 are likely to operate together. 

83 Both instabilities occur more slowly than thermal equilibration of the atmosphere and surface. This 

84 separation of timescales allows us to solve for the fast processes that set the surface temperature 

85 (in §2), and then separately address each of the two slower processes which may cause atmospheric 

86 pressure to change (in §3-§4). 

87 Day-night color temperature contrast is among the 'easiest' parameters to be measured for a 

, but there is currently no theory for ATg on planets with 



89 observable surfaces. One motivation for this paper is to contribute to this emerging theory. We use 

90 a simplified approach which complements more sophisticated exoplanet Global Circulation Models 



j[GCMs) (Joshi et al. 


1997 


Joshi 


2003 


Merlis &: Schneider 


2010 


Edson et al. 


2011 


Wordsworth 


|et al.| 


|2011 


Pierrehumbert 


2011) 


. §2 sets out the energy balance model that is used for both 



93 instabilities, and §3 explains the enhanced substellar weathering instability including our choice of 

94 weathering parameterization. §4 explains the substellar dissolution feedback (considering only the 

95 1:1 spin-orbit resonance). We find that SDF doesn't work in most cases, so readers motivated by 

96 short-term detectability can safely omit §4 and move to §5. §5.1 discusses relevant solar system 

97 data, including the possibility that Mars underwent a form of ESWI. §5.2 discusses applicability to 

98 exoplanets, and §6 summarizes results. 



transiting exoplanet (Cowan &: Agol|2011 



2. Idealized energy balance model 



Consider a planet in synchronous rotation on which surface liquid water is stable, with an at- 
mospheric temperature that decreases with height along the dry adiabat. Slow rotation weakens 
the Coriolis effect, allowing the atmospheric circulation to all but eliminate horizontal gradients 
in atmospheric temperature at the top of the boundary layer, T^. This is the weak temperature 



gradient approximation often made for Earth's tropics (e.g. Merlis & Schneider (2010)). Figure 
[T] shows the setup for our idealized energy balance model. The surface temperature Ts{ip) at an 
angular separation ifj from the substellar point is set by the local surface energy balance:- 



SWs{il^) - LW\{ilj) + LWi -/3(T,(V) - Ta) = 



(1) 



where SWs{'4') is starlight absorbed by the surface, LW[{iIj) = r]aTs{tlj)^ (where a is the Stefan- 
Boltzmann constant and ?] ~ 1.0 is the emissivity at thermal wavelengths) is the surface thermal 
radiation, LW^ is backradiation from the atmosphere, /3 is a turbulent heat transfer coefficient (/3 = 
kTF P, where p is the near-surface atmospheric density divided by Earth's sea- level atmospheric 
density, and kxF is a turbulent flux proportionality constant), and Tq is the temperature of the 
atmosphere at the top of the boundary layer. (Equatorial superrotating jets can cause the hottest 



113 ipoint on the surface to be downwind from the substellar point (Knutson et al. 2009; Mitchell & 



114 Vallis 2010; Liu & Schneider 2011)). The shortwave flux SWs{ip) = L*(l — a)cos(V') corresponds 



- 5 - 



115 to stellar flux L*, attenuated by surface albedo a. There is negligible transport of heat below the 

116 surface: we assume seas are not globally interconnected (or are shallow, or are deeply buried, or do 

117 not exist) and energy flux from the interior is small. 

118 LW\^ ( = \ ^^LW\i(i\)) sin dV' — OLR, ignoring turbulent fluxes) is longwave flux from the at- 

119 mosphere to the surface. OLR (Outgoing Longwave Radiation, longwave energy exiting the top of 

120 the atmosphere) is given by interpolation in a look-up table. To build this look-up table, we slightly 

121 modified R.T. Pierrehumbert's scripts at http://geosci.uchicago.edu/~rtpl/Principles 

122 PlanetaryClimate/, particularly PureC02LR.py. The look-up table gives OLR(P\,Ts) for a pure 

123 noncondensing CO2 atmosphere on the dry adiabat, with the temperature at the bottom of the 

124 adiabat equal to the energy- weighted average Tg, and with Earth gravity (9.8 m/s^). Our noncon- 

125 densing assumption introduces large errors for Tg < 175K, so we assume that the top-of-atmosphere 

126 effective emissivity OLR/LW'[ at these low temperatures is the same as at Tg = 175K. 

127 To investigate atmospheres not made of pure CO2, we introduce an opacity ratio or relative 

128 radiative efficiency A, which is the ratio of the radiative efficiency of the atmosphere of interest 

129 to that of pure noncondensing 002- A is a simplification of the complicated behavior of real gas 

130 mixtures ( Pierrehumbert |2010 ). A can be greater than 1 if the atmosphere contains radiatively 



131 very efficient gases (chloroflourocarbons, CH4, NH3, or the "terraforming gases"; Marinova et al 



132 (2005)). We then query the look-up table using Pa = AP. Smaller values of A have a weaker 



133 greenhouse effect (increased OLR). 

134 Rayleigh scattering is relatively unimportant for planets orbiting M-dwarfs. Starlight is con- 

135 centrated at red wavelengths where Rayleigh scattering is ineffective (falling off as A~^, where A is 

136 wavelength). The optical depth to Rayleigh scattering of 1 bar of Earth air is 0.16 for light from 

137 the Sun, but only 0.02 for light from the Super-Earth hosting M3 dwarf Gliese 581 (approximating 

138 both stars as blackbodies) . In addition, absorption of starlight by the atmosphere is much stronger 

139 in the NIR than the visible, and so is more effective at compensating for Rayleigh scattering as 



140 star temperature decreases (Pierrehumbert 2010). We neglect Rayleigh scattering and absorption 

141 of starlight by the atmosphere. 

142 The horizontally uniform atmospheric boundary layer temperature, Tq, is set by the total 

143 energy balance of the atmosphere. 



If [LW-[ {7P) + /3{Ts{iP)-Ta)] sin iP dip -OLR-LWi=0 (2) 
^ Jo 

144 where the integral gives the average flux from the surface. This reduces to = 1/2 J^Tg^ip) sin tp dip 

145 because of our particular choice of LW],. In effect, we assume that the boundary layer only interacts 

146 with the ground through turbulent fluxes. 

147 For a given Pa, we iterate to find and Ts{ip). ip resolution is 5 degrees. The initial 

148 condition has the surface in radiative equilibrium, and the atmosphere in equilibrium with this 
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149 surface temperature distribution. Convergence tolerance is ~ 2 x 10 



150 Throughout the paper, we assume kxp = Cp{Ta)CDU, where CpiTa) is the temperature- 

151 dependent specific heat capacity of CO2 (~850 J/kg/K at 300K), C/j is a drag coefficient, and U is 

152 characteristic near-surface wind speed. Cd 



fr2 



where kyx = 0.4 is von Karman's constant. 



153 zi = 10m is a reference altitude, and zq is the surface roughness length (10~^ m, which is bracketed 

154 by the measured values for sand, snow, and smooth mud fiats; Pielke (2002)). For our reference 

155 U = IQ m/s, this gives kxp P = 12-3 W/m^/K. We assume a Prandtl number near unity. 

156 reports sensitivity tests using different values of U. 



55.2.1 



157 We accept the following inconsistencies to reduce the complexity of the model:- (1) Radiative 

158 disequilibrium drives Ts to a higher value than Tq at the surface, and turbulence can never com- 

159 pletely remove this difference. Therefore, setting the temperature at the bottom of the atmosphere 

160 to the surface temperature will lead to an overestimate of LWX at low P. (2) The expression for kxp 

161 is appropriate for a neutrally stable atmospheric surface layer, but turbulent mixing is inhibited on 

162 the nightside by a thermal inversion ( Merlis Schneider|2010 ) . This will tend to make the coupling 

163 between nightside atmosphere and nightside surface too strong. (3) We assume an all-troposphere 



atmosphere with horizontally uniform temperature. Merlis & Schneider (2010) find temperatures 
are nearly horizontally uniform for Earthlike surface pressure and for levels in the atmosphere at 
pressures less than half the surface pressure. Atmospheric temperature is approximately horizon- 
tally uniform when the transit time for a parcel of gas across the nightside, Tadvect = 7^1 is short 
compared to the nightside radiative timescale, Trad 



p Cp 



169 planet radius (1 Earth radius), Uh is high-altitude wind speeci(^ 



Showman et al. 


2010 


). Here, a is 


i (~ 30 m/s: M 


erlis & Schneider 



170 (2010)), e is an greenhouse parameter corresponding to the fraction of the emitted radiation that 

171 is not absorbed by the upper atmosphere and escapes to space, and T = 250K is the atmospheric 

172 temperature, the radiative equilbrium temperature on the darkside being zero. Picking e = 0.5, 

173 this gives Tadvect ~ 2 days and T^ad ~ 50 days. (4) The treatment of Tq is crude. (5) We assume 

174 the atmosphere is transparent to stellar radiation, which is a crude approximation under M-dwarf 

175 (or cloudy) skies. (6) We neglect condensation within the atmosphere. 

176 Representative temperature plots are shown in Figure § for A = 0.1 and = 900 W/m^. 

177 At low pressures, nightside temperatures are close to absolute zero, and substellar temperature is 

178 close to radiative equilibrium. Increasing pressure cools V < 60°, and warms ip > 70°. This is 

179 beause the atmosphere is warmer than the surface on the nightside, and cooler than the surface 

180 on the dayside. Therefore, the increase in P(oc (3) increases the f3{Ts — Ta) term, which warms the 

181 nightside, but cools the dayside. For positive A, LW^ will increase with P and warm the entire 

182 planet. However, for the relatively small value of A shown here, the substellar point still undergoes 

183 net cooling with increasing pressure. This cooling with increasing P is what makes the ESWI and 

184 SDF possible. When the surface becomes nearly isothermal, as for the "10-bar" curve in Figure [2j 

185 the entire surface warms with increasing pressure, and the ESWI and SDF cannot occur. 



For A > 0, Tg must increase with P. Even if there is no greenhouse effect, the homogenization 
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187 of the atmosphere wih warm the planet on average because of the nonhnear dependence of T on 

188 energy input (Edson et ah 2011). For small optical depth, nightside Tg oc ln(P). The fractional 

189 area of the planet over which liquid water is stable is 27% at radiative equilibrium, decreasing with 

190 pressure and vanishing at ~0.7 bars. As P increases the greenhouse effect further, liquid stability 

191 reappears at ~2.4 bars, rapidly becoming global. 



192 
193 



3. Climate destabilization mechanism #1: 
Enhanced substellar weathering instability (ESWI) 

3.1. Weathering parameterisation 



196 
197 



The Berner Sz Kothavala (2001) weathering relation, which is specific to CO2 weathering of Ca-Mg 



silicate rocks, states 



Wo 



P 



0.5 



exp [kACT (TsW - To)] [1 + knuN (TsW - To)] 



0.65 



(3) 



direct T dependence 



hydrology dependence 



where is local weathering rate, Wq is a reference weathering rate, P is atmospheric pressure, Po 
is a reference pressure. To = 273K is a reference temperature, kACT = 0.09 is an activation energy 
coefficient, and knuN is a temperature- runoff coefficient fit to Earth GCMs. Equation 3 is widely 
used, but uncertain and controversial (§5.2.1). In our model, the strong temperature dependence 
leads to a strong concentration of weathering near the substellar point. For example, for a 1-bar 
atmosphere with A = 0.1 (shown in Figure [2]), 93% of the weathering occurs in 10% of the planet's 
area. The planet-averaged weathering rate is 



Wt{P) 



1 

47r 



W^A^ dV' 



(4) 



Climate is in equilibrium (^=0) when planet-integrated weathering of greenhouse gases, Wt, is 
equal to net supply Vn by other processes. The climate equilibrium is stable if dWt/dP > - 



in this case, carbonate-silicate feedback enables long-term climate stability (Walker et al. 1981). 



Climate destabilization occurs when dWt/dP < 0. The weathering feedback that underpins the 



Circumstellar Habitable Zone concept (Kasting et al. 1993) changes sign, and acts to destabilize 
these climates. 
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3.2. ESWI results 



212 Figure [3] shows weathering rates corresponding to the temperatures in Figure [2} The hues show 

213 equihbria between weathering consumption of greenhouse gases and net supply by other processes 

214 (Figure [3]). The units of weathering are normahzed to the weathering rate at P = 3 mbar. Along the 

215 curves, planet-integrated weathering of greenhouse gases, Wt, is equal to net supply Vn by other 

216 processes. These other processes can include volcanic degassing, metamorphic release, biology, 

217 sediment dissolution, and loss to space. The shape of this curve is set by competition between 

218 three effects:- (1) the greenhouse effect (A), (2) advective heat transport (/3,/ctf) which enlarges 

219 and maintains unstable regions, and (3) stellar flux (L^,), which shrinks the unstable region. The 

220 curve has two stable branches and an unstable branch. The slope of the low-pressure stable branch 

221 is set by the (P/Pq) term in Equation 5 - for small P, t oc ln(P) and /3 oc P are small and the 

222 atmosphere has little effect on Tg . On the intermediate pressure unstable branch - dashed in Figure 

223 [2]- the atmosphere is important to energy balance but /3 outcompetes r. The homogenizing effect of 

224 f3 cools the substellar region, and planet-integrated weathering decreases with increasing pressure. 

225 In the high-pressure stable branch, the planet is close to isothermal. Further increases in (3 have 

226 little effect on ATg, but r warms the whole planet and now is able to increase Wt- Catastrophic 

227 jumps in the pressure caused by small changes in the supply Vn occur at ~0.05 bar (for increasing 

228 volcanic activity) and ~1 bars (for decreasing volcanic activity). The jumps correspond to a >10^- 

229 fold increase in pressure, or a > 10^-fold decrease in pressure, respectively. The existence and 

230 location of these bifurcations are sensitive to small changes in the coefficients of (3). The timescale 

231 for the climate regime jump is set by the rate of weathering and/or rate of volcanism on each 

232 specific planet. For example, Earth today supplies ~15 mbar CO2 in 10^ yr (atmosphere + ocean: 



233 linearizing, ~2 x 10 yr to build up 1 bar CO2), but an lo-like rate of resurfacing (Rathbun et al. 

234 2004) with the same magmatic volatile content would build up 1 bar in ~ O{10^-^) yr. A natural 



weathering-rate experiment occurred on Earth 0.054 Gya, with very rapid release of CO2 from an 
unknown source. The warmed climate required At ~ O(IO^) yr (Murphy et al. (2010), using '^He 



237 accumulation dating) to draw down 0.9 mbar of excess CO2 (Zeebe et al. 2009). 



238 Figure |4] shows habitable-zone climate regimes as a function of equilibrium pressure and stellar 

239 flux, which can change (Figure |4^) due to stellar evolution, tidal migration (Jackson et al. 2010), 

240 or close encounters with other planets and small bodies (Morbidelli et al. 2007). Climate stability 



241 depends strongly on A, so we show this for three values of A - an almost radiatively inert gas (Figure 

242 ^p), an intermediate A = 0.1 case (Figure [4|;), and a strong opacity that only just allows the ESWI 

243 (Figure [4]i, for A = 1.0; higher values are stable against ESWI). The thick black line labelled with 

244 zeros corresponds to marginally stable climate equilibrium. Increasing L=k widens the range of P 

245 that sits within the low-pressure stable branch. That is because higher produces higher absolute 

246 temperatures, and higher absolute temperatures favor radiative exchange between atmosphere and 

247 surface (oc 4(Ts — Ta)"^ for small (T, — To) ) versus turbulent exchange (oc (T^ — T^)^). Therefore, a 

248 given equilibrium value of P on the low pressure branch is more stable at low than high :- for a 

249 small increase in P the radiative warming will be less counteracted by the cooling of the substellar 



- 9 - 



point. Increasing pushes the the high pressure branch to higher values. Other instabihties 
that contain the habitable zone are shown by thin lines. The dash-dotted line corresponds to 
pressures in excess of the nightside CO2 saturation vapor pressure. CO2 atmospheres to the left of 
this line condense on fast, dynamical timescales. Increasing A couples the horizontally-isothermal 
atmosphere more strongly to the nightside surface, and causes the CO2 collapse threshhold to move 
to the left. The thin solid line corresponds to mean surface temperatures above 40° C, which is the 
lower edge of the moist runaway greenhouse zone ( Kastingl|1988 Pierrehumbert 1995) - our model 



does not include fluxes of latent heat or the greenhouse effect of water vapor, so the position of 
this line is notional. The maximum pressure and minimum pressure of the bifurcation loop (Figure 
p3h are shown by thick gray lines in figure m In some cases, the lower end-of-transition pressures 
are so low that water boils (thin dashed line). This would initiate the boiling of a global ocean. If 
boiling continued, the eventual fate could be a steam ocean, or restriction of liquid water stability 
to a thin belt near the terminator. The upper end-of-transition pressures are often >10 bars, with 
a nearly isothermal Tg. Such planets would have weak and perhaps undetectable phase curves in 



264 the absorbing bands (Selsis et al. 2011). 



For the almost radiatively inert gas (Figure |4j3), geologically stable equilibria usually have 
P <0.01 bars or P >0.3 bar. Intermediate pressures cannot be stable over geological time. The 
overall pattern is similar for A = 0.3 (Figure |4|:) . Increasing A always shrinks the domain of the 
unstable branch. For much higher A the climate is stable everywhere. We show the last gasp of 
the instability in Figure |4]i. Higher radiative efficiency means that for a small change in P, when 
AT5 is significant, radiative heat transfer overcomes advective heat transfer, the substellar patch 
warms, and overall planet temperature increases. 

Figure [s] summarizes the effects of ESWI in a stability phase diagram (against axes of gas 
radiative efficiency and incident stellar radiation). For A > ~1, the climate is stable to ESWI. For 
A < ~1, ESWI is possible but the pressure jumps caused by ESWI do not always have a catastrophic 
effect. Higher L^, warms the climate towards the runaway moist greenhouse threshhold, and upward 
jumps in pressure for L^, > ~2000 W m~^ may initiate the moist runaway greenhouse (points to 
the right of the vertical dashed line in Figure 5). Atmospheric collapse to ~ 1 mbar (well below the 
triple point of water) only occurs below A < 0.4. Increasing L^, increases both the bottom and the 
top pressure for instability, implying that tidal migration towards the star should be destabilizing 
for thick atmospheres but stabilizing for thin atmospheres. 



4. Climate destabilization mechanism #2: Substellar dissolution feedback (SDF) 

Water ice and basalt are the most common planetary surface materials in the Solar System, and 
are expected to be common elsewhere. When these melt (around 273K for ice and 1300K for basalt), 
atmosphere can dissolve into the melt. Counterintuitively, a decrease in P and in average surface 



285 temperature T, can favor melting if AT^ is large, as pointed out for Mars by Richardson & Mischna 



286 (2005). For a synchronously-rotating planet entirely coated in condensed material (ice, rock, or 
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287 carbon-rich ceramic), surface liquid will first appear near the substellar point. Atmosphere will 

288 dissolve into this warm little pond, approaching Henry's-law equilibrium: 



Ppond — j-j ( D 

^Earth 



pond 2 ^'^^ ^max , 



P^mpi kHiT")exp 



C 



T, 



pond 



1 



(5) 



global— equivalent liquid depth 



mass of gas per unit liquid volume 



289 where Ppond (h^ bars) is the equivalent atmospheric pressure of gases dissolved in the ocean, g is 

290 surface gravity, PEarth = 1-01 x 10^ Pa is a normalization constant, Dpond is pond depth, pi is 

291 liquid density, il^max is the angular radius of the pond, P is the surface pressure, d is a dissolution 

292 exponent (~0.5 for water in silicate liquids and ~1 for gases in water), m is the molecular mass 

293 of the atmosphere, kn and C are Henry's-law coefficients, Tpond is the pond temperature, and T° 

294 is a reference temperature. Here, the first term in brackets is the depth of the pond in global- 

295 equivalent meters, and the second term in brackets is the mass of gas per unit volume of pond. 

296 We have neglected the distinction between fugacity and partial pressure. The pond is assumed to 

297 be well-mixed so that Tpo^^ = j^"^""^ Tg^^ slnil: dil) . This relation assumes a uniform heat 

298 transfer coefficient between the surface and the pond. 

299 Instability occurs when a decrease (increase) in surface pressure results in an uptake (release) 

300 of gases from the pond that exceeds that which is cons istent with that change in pressure, i.e. 

301 



^ Ppond 



dP 



< — 1. In this case the feedback has infinite gain (Roe 2009). Pond growth rate is limited 



302 by the balance between insolation and the latent heat of melting (O(IO^) yr for melting of a 1km- 

303 thick ice sheet. Earthlike insolation, a = 0.6, and 10% of sunlight going to melting), or by thermal 

304 diffusion of heat into a stratified pond. For tectonic activity levels that are not too high, the SDF 

305 is much faster than the ESWI, which is limited by volcanic degassing rates and tectonic resurfacing 



306 rates (0.05 mm/yr on Earth (Leeder 1999); 16 mm/yr on lo (Rathbun et al. 2004)). The SDF 

307 stops when insolation is insufficient to allow further pond growth. What happens after the SDF 

308 has operated will depend on the sign of the carbonate-silicate feedback at the new, modifed P. If 

309 dWt/dP > 0, the normal carbonate-silicate feedback will rejuvenate the atmosphere on a volcanic 

310 degassing timescale, freezing the ocean. If dWt/dP < 0, pressure decreases further. 



Gases that react chemically with seawater (such as SO2 and CO2; Zeebe & Wolf-Gladrow 



312 (2001)) can have Ppond niuch greater than that given by Henry's law. For example, total Dissolved 

313 Inorganic Carbon (DIG, oc Ppond) is buffered against changes in P by carbonate chemistry, and 
3w Ppond changes much more slowly than Henry's law. Ppond oc P^'^ for the modern Earth ocean, as 



315 opposed to Ppond oc P for Henry's law iZeehe k Wolf-Gladrow 2001 Goodwin et al. 2009^. For 

316 a decrease in P that increases ipmax, this buffering favors the tendency of increasing pond volume 

317 to draw down more atmosphere, against the Henry's Law decrease in atmospheric concentration 

318 per unit pond. Garbonate buffering is less important for P >~1 bar, because at the correspond- 

319 ingly low pH the DIG partitions almost entirely into GO2. For CO2, we use R. Zeebe's scripts 

320 (http : //www. soest .hawaii . edu/oceanography/f aculty/zeebe_f iles/C02_Systein_in_Seawater 
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/csys.html) to find the additional DIG held in the ocean as HCO3 and CO3 . We use fixed alka- 
linity, 2400 //mol/kg (similar to the present Earth ocean; Zeebe &: Wolf-Gladrow ( 2001[ )). Figurejo] 



323 shows the results, ipmax is set by Equation 1 through Ts{ip). Ocean circulation adds additional heat 

324 transport terms to Equation 1, which we ignore. We also do not consider buffering by dissolution 

325 and precipitation of carbonates or salts (such as sulfates). 

326 As an example of SDF, consider the partitioning of an initial 1-bar total inventory of CO2 

327 (blue-green contour labelled '0' in Figure [g]) with increasing L^,. Initially, the planet's surface is 

328 below freezing. As stellar flux is increased to 800 W/m^, a small sea forms and dissolves some of 

329 the CO2. The system is within the area where SDF is a positive feedback on small changes in P 

330 (outermost black contour in Figure [g]) , which accelerates sea growth. A small further increase in 

331 leads to runaway SDF (innermost black contour in Figure [g]) , and ocean area quickly grows with 

332 from ~5% to ~15% of planet surface area with no change in L*. After this change the atmospheric 

333 inventory is reduced to 1/4 bar with the remaining 3/4 bar dissolved in the ocean. Further pond 

334 growth requires further increase in L*. Increasing Tpond decreases solubility and shallows the slope 

335 of increasing volatile inventory stored in the ocean. By L^, = 3800 W/m^ (the highest considered), 

336 the ocean covers almost the entire lightside hemisphere, and stores ~ Ye of the initial CO2 inventory. 

337 Further small increases in L^, cause a global ocean. 

338 Three main effects control this system:- (1) Cosine fall-off of starlight weakens the ability of 

339 decreasing pressure to increase ocean area beyond a relatively small Vmax- Following a line of 

340 decreasing pressure, the dashed red lines (fractional ocean area) become more widely spaced with 

341 decreasing pressure. Cosine falloff of stellar radiation is responsible. This restricts the scope of 



342 the instability, which tends to lead to Eyeball states ( Pierrehumbert 2011). We do not find any 

343 cases where the SDF can turn a dry planet into an ocean-covered planet or vice versa. (2) Ocean 

344 instability disappears above ~5 bars, when the surface is nearly isothermal. Because decreases in 

345 pressure always decrease Tg (except when there is an antigreenhouse), decreases in pressure can only 

346 cause oceans to freeze over. Below ~5 bars, AT^ is not negligible. If a decrease in pressure allows 

347 the substellar temperature to rise above freezing, an ocean can form. (3) Rectification of starlight 

348 by the terminator, and of ocean area by the melting point, divides the phase diagram into 'no 

349 ocean', 'substellar ocean' and 'global ocean' zones. On the nightside T, is constant, so ocean area 

350 jumps from 50% to 100%. Notice the change in sign of dependence below and above the dashed 

351 red line corresponding to 50% ocean area. In the substellar ocean zone, increasing increases 

352 ocean area and the ocean inventory increases. However, in the global ocean zone, increasing 

353 cannot increase ocean area. The decrease in gas solubility with increasing temperature dominates, 

354 and the ocean inventory decreases. 

355 Figure |6] does show unstable regions in parameter space for which substellar dissolution is a 

356 positive feedback on changes in P (solid black lines). These always correspond to small oceans 

357 (<10% of planet surface area). But the gain of the feedback is small, and runaways cannot occur 

358 unless ocean depth >10km. Factor-of-3 decreases (or increases) in atmospheric pressure can occur 

359 with no change in the total (atmosphere + ocean) inventory of volatile substance. The SDF is a 
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minor positive feedback (or a negative feedback) for most gases which dissolve simply according 
to Henry's Law. The buffering effect is not sufficient to allow ocean area to give a strong positive 
feedback. Finally, A = 0.03 is unrealistically low for all-C02 atmospheres. Setting A = 1 would 
shut down the instability. We conclude that the CO2 SDF is unlikely to be important for planets 



jn synchronous rotation, and is important in fewer cases than is the ice-albedo feedback (Roe & 
Baker||2010l [Pierrehumbert et al. |[2011[ ). 



Discussion 



5.1. Solar system climate stability - is Mars a solar system example of ESWI? 

Although ESWI is most relevant for planets in synchronous rotation, it can work for any planet 
with sufficiently high ATg. Therefore, in order to test ESWI against available data, we compare the 
requirements for ESWI to Solar System data (Table 1). After correcting for the distorting effects of 
life, all of the Solar System's non-giant atmospheres are overwhelmingly one gas. Except for Earth, 
the principal gas is also the main greenhouse gas. Venus' atmospheric composition is not controlled 



py the abundance of surface liquid (nor in solid-state equilibrium with surface minerals; Hashimoto 
& Abe (2005); Tremain &; Bullock ( |in press )), and Triton's atmosphere is too thin to stabilize 
liquid nitrogen. Climate regulation on Titan is not well understood. Currently, the greenhouse 



effect of CH4 outcompetes the antigreenhouse effect of the organic haze layer ( McKay et al.|[l99T ) 



377 The production rate of the organic haze layer depends on [CH4] (McKay et al. 1991; Lorenz et al. 



1997). ESWI is not currently possible on Titan because ATs is too small, 2.5-3.5K (Jennings et al. 



2009). Therefore, out of 5 nearby worlds with atmospheres and surfaces, only Mars is a candidate 
for ESWI (§5.1). Solar System data suggest that the conditions for the ESWI are quite restrictive, 
and that most planets will not be susceptible to the ESWI. 

Mars may have passed through ESWI in the past. The main climate-controlling greenhouse 
gas (CO2) can dissolve in liquid water and be sequestered as carbonate minerals. Mars has low P 



384 (95% CO2). It has widespread deposits of carbonate minerals (Bandfield et al. 2003 Wray et al. 



2011 ), but little or no surface liquid water. It currently sits at the gas- liquid sublimation boundary 



386 (Kahn 1985), with ATs ~100K , and GCMs show that AT^ as Pf (Richardson k Mischna 2005). 



These geologic and climatic observations are all consistent with a past rapid transition via ESWI 
from an early thicker atmosphere to the current state, as follows. Increasing solar luminosity could 
have permitted transient liquid water, allowing carbonate formation. The corresponding drawdown 
in P would increase noontime temperature, allowing a further increase in liquid water availability 
and the rate of carbonate formation. The runaway would slow as P approached the triple-point 
buffer: loss of liquid water stability may have throttled weathering and buffered the climate near 



393 the triple point (Kahn 1985 Halevy et al. 2007). However, ice is not ubiquitous on the surface 



394 and can migrate away from warm spots, so unusual orbital conditions are neccessary for melting 
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( Kite et al.||2011a|b ). Alongside carbonate formation, atmospheric escape, polar cold trapping and 
volcanic degassing are the four main processes affecting P on Mars over the last 3 Ga. However, 



397 recent volcanism has been sluggish and probably C02-poor (Hirschmann & Withers 2008; Werner 



2009 Stanley et al. in press), and present-day atmospheric escape appears slow (Barabash et al. 



2007). Polar cold traps hold ~1 Mars atmosphere of CO2 as ice today (Phillips et al. 2011), but 



this trapping should be reversable at high obliquity. Therefore, it is quite possible that substellar 
(= noontime) carbonate formation has been the dominant flux affecting the secular evolution of 
Mars' atmosphere since 3.0 Gya. This hypothesis deserves further investigation. 
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5.2. Application to exoplanets 

5.2.1. How general is our feedback? 



405 ESWI requires: 



- Strong temperature dependence of the weathering drawdown of greenhouse gases. Decreasing 
^ACT to 0.03 (from the nominal 0.09) eliminates the ESWI except for radiatively inert atmospheric 
compositions. On the other hand, increasing k^cT to 0.27 causes a large unstable region even for A 
= 1.0, with at least a 1 dex range of atmospheric pressure unstable to ESWI for all habitable-zone 
luminosities. 

Earth data on the value of k^cT are consistent with k^cT ~ 0.1. However, deep-time calibra- 
tion of weathering-temperature relations such as Equation 3 is difficult. There is only one planet 
(Earth) to use as laboratory, with constantly drifting boundary conditions, and rather few natural 
experiments. Chemical weathering rates of silicate minerals in the lab are definitely temperature- 



dependent (White & Brantley 1995), but erosion-rate dependence is also important at the scale 



of river catchments (West et al. 2005). Confirming temperature dependence on geological scales 



is difficult, in part because today's weathering rate contains echoes of glacial-interglacial cycles 



( Vance et al. 2009 ) . Regression of present-day river loads on watershed climatology by West et al 

^®'^0s/^^'^0s data suggest continental 



(2005) suggests an e-folding temperature of 8.5(+5.5/-2.9)K. 
weathering rates increased 4-8x in <ClO^ yr during a Jurassic hyperthermal (AT < lOK) 0.183 Gya 
(Cohen et al. 2004), implying an e-folding temperature <(5 ~ 7)K. Analysis of the apparent time 



dependence of weathering rate gives support to a hydrological control on weathering rates ( Maher 



2010 


), but on c 


Schneider 


2008 



DW 
Dt 



; rates ( 


Maher 


0' Gorman & 



gOlO), but on a planetary scale precipitation always increases when Tg increases (O'Gorman & 
■ ■ ~— dw I dw dR -.1, an ^ n 



Overall, Equation 3 is consistent with deep-time, present-epoch, and laboratory estimates for 
Earth. Though Equation 3 is used in this paper as a general rule for Earth-like planets, the 
weathering-temperature relation is shaped by biological innovations. For example, the symbiosis 
between vascular plants and root fungi (arbuscular mycorrhizae) acidifies soil, profoundly accel- 
erates weathering, and may be unique to Earth ( [Taylor et al. 2009). All geologically important 
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430 surface weathering reactions require a liquid phase ( [White fc Brantley|[T995| ) . We assume weather- 

431 ing reactions do continue below the freezing point (due to microclimates, or monolayers of water). 

432 However, the results shown in Figure |4] and [5] did not change qualitatively when we set W = for 

433 T < Tmelt- 
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Strong temperature dependence could break down in several ways. For example, erosion is 
needed to expose fresh mineral surfaces for weathering. On a tectonically quiescent planet (and 
for tectonically quiescent regions of an active planet) weathering may be limited by supply of fresh 
surfaces, with weathering going to completion on all exposed silicate minerals. On tectonically very 
active planets with little or no land, hydrothermal alteration and rapid seafloor spreading maintain 
low greenhouse gas levels with little or no temperature dependence (Sleep Sz Zahnle [2001 ). A 
planet with weathering rates that are limited by the availability of liquid water (like Mars) can 
have effectively temperature-dependent weathering, but only if the atmospheric pressure is well 
above that fluid's triple point. 

- Large ATg. The high-AT^ requirement cannot be met if the atmosphere is thick. A deep global 
ocean circulation behaves like a thick atmosphere - Earth abyssal temperatures vary <2K from 
tropics to poles ( Schlitze"r||2000 ). Therefore, the climate destabilization mechanism cannot operate 
on a planet with a deep global ocean circulation - at least some land is necessary. Evaporation 
will dry out land at the substellar point if Tg is high, so weathering activity may be concentrated 



ipt cooler il> in this case. ATg varies little, or even increases, with rotation frequency (Merlis & 



Schneider 2010 Edson et al. 2011). Therefore, ESWI could work for rapidly-rotating planets such 



as Mars (§5.1). However, the isothermal approximation does not apply when the Coriolis force 
prevents fast equator-to-pole winds. For rapid rotators, Tg is a function of latitude and longitude, 
and our idealized energy balance model is not appropriate. 

- Small A. Strong greenhouse gases have high A, which suppresses ESWI. On the other hand, 
A can be negative if there is an antigreenhouse effect (A < 0). M-dwarfs later than M4, with 
fully convective interiors, seem to remain active with high UV fluxes for much longer than do 



Sunlike stars. High UV fluxes broadly favor CH4 accumulation (Segura et al. 2005) and perhaps 



antigreenhouse haze effects. When A < 0, ESWI will apply for all P and L^. 

- Surface- atmosphere coupling. This consists of radiative and turbulent coupling. Turbulent cou- 
pling requires a nonzero near-surface wind speed, and that the global near-surface atmosphere is 
not stably stratified. A sensitivity test with a 10-fold reduction in U showed that the pressure 
range unstable to ESWI moves to ~ 10 x higher pressure. The range of A subject to ESWI was 
significantly reduced. We assume U is not a function of P, but simulations show that U increases 
as P decreases (Melinda A. Kahre, via email). This would strengthen the instability. 

SDF has very similar requirements to ESWI, but the constraint on rotation rate is stricter. Large- 
amplitude libration or nonsynchronous rotation would prevent the development of a deep pond 
around the substellar point. (A low-latitude liquid belt can be imagined, but the idealized EBM 
of §2 is not appropriate to that case). Kepler data show that only a small proportion of close- 
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468 in small-radius exoplanets in multi- planet systems are in mean-motion resonance (Lissauer et al. 



2011), but that most are close to resonance and could maintain nonzero eccentricity. This would 



allow for significant nonsynchronous rotation if the planet's spin rate adjusts to keep the substellar 
point aligned with the star during periapse passage (pseudo-synchronous rotation). In addition, 
SDF requires very soluble gases: by contrast, the N2 content of even a 100km deep ocean at 298K 
is only ~0.2 bar per bar of atmospheric N2. In this paper we assume that volatiles are excluded 
from the ice when the ice freezes. If clathrate phases form, they could absorb volatiles and make 
the SDF irreversible. 



5.2.2. Climate evolution into the unstable region 



477 Planets could undergo ESWI early in their history if they form in the unstable region. In addition, 

478 many common geodynamic and astronomical processes can shift the equilibrium between Wt and 

479 Vn (Figure |4^), causing a secular drift of the equilibrium across the phase space {L^,, A, Vt} (Figures 

480 [3] and [5]). This drift can take a planet from a stable equilibrium to an unstable equilibrium via a 

481 saddle-node bifurcation (Strogatz et al. 1994). 



• Dynamics and stellar evolution: Theory predicts that the secular increase in solar flux should 
have gradually shifted the position of Earth's climate equilibrium. This is consistent with the 
sedimentary record of the last 2.5 x 10^ yr ( Grotzinger &: Kasting|[T993 Grotzinger fc James||2000 



Ridgwell & Zeebe 2005 Kah & Riding 2007). For Kepler field (Sunhke, rapidly evolving) stars, 
this could also occur and potentially cause the ESWI for planets with initially thick atmospheres 
(Figure 4a). For M stars, main-sequence insolation changes little over the lifetime of the Universe. 

• Atmospheric evolution: A can change as atmospheric composition evolves. For example, the rise 
in atmospheric oxygen following the emergence of oxygenic photosynthesisers probably oxidized 



490 atmospheric CH4 and may have caused a catastrophic decline in A (Kopp et al. 2005 Domagal- 



Goldman et al. 


2008 


. Even 


broadening ( 


Li et al. 


2009 


)• 



where atmospheric erosion exceeds geological degassing. For these planets, Vn is negative. Strong 
stellar winds and high XUV flux are observed for many M stars. Removal of atmosphere by strong 
stellar winds (Mura et al. 2011) or, for smaller planets, high XUV flux ( |Tian 2009) could trigger 
an instability for a planet orbiting an M-star, by reducing P (Figure 4a). 



Tectonics and volcanism: Volcanic activity decays with radioactivity ( Kite et al.|2009 Sleep|2000 



2007 Stevenson 2003), so in the absence of tidal heating the equilibrium pressure will gradually 
fall (on a stable branch where Wt increases with P) (Figure 4a). Superimposed on this decline 
are pulses in volcanic activity due to mantle plumes, and perhaps planetwide volcanic overturns 
as seem to have occurred on Venus. This will cause spikes in equilibrium pressure. The rate of 



502 resurfacing is very sensitive to mantle composition, tidal heating, and tectonic style (Kite et al 



g009l [Valencia U'Connen|[2009l |Korenaga|[20Tol [Behounkova et al.|[2010l [20Tl| [van Summeren et 
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ar]|2011 )■ Shutdown of volcanism (such that Vn < 0) extinguishes the possibihty of a stable chmate 



equihbrium; P will fall monotonically. ESWI can accelerate this decay, and Mars may be an example 
of this (§5.1). Mountain range uplift exposes fresh rock and may provide a O(IO^) yr increase in 
weathering rate that cools the planet (as arguably and controversially may have occurred for Tibet, 



Earth: Garzione (2008)). This may trigger ESWI by lowering pressure. More speculatively, drift 
of continents can cause very large atmospheric pressure fluctuations if greenhouse-gas drawdown 
occurs mainly on land and is strongly focussed in a high-weathering patch near tp = 0. A continent 
drifting over the patch will increase planet-averaged weatherability, and pressure will go down. 



6. Summary and conclusions 
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Nearby M-dwarfs are targeted by several planet searches:- MEarth (Charbonneau et al. 2009); the 



VLT4-UVES M-dwarf planet search (Zechmeister et al. 2009); the VLT+CRIRES M-dwarf planet 



search (Bean et al. 2010); HARPS (Forveille et al. 2011); M2K (Apps et al. 2010); and proposed 



space missions TESS ( |Deming|[2009l ) ; ELEKTRA; PLATO; and ExoplanetSat ( [Smith et al.][2010| ). 
These searches are driven in part by the hope that planets orbiting M-dwarves can maintain surface 
liquid water and be habitable. Maintaining surface liquid water over geological time involves 
equilibrium between greenhouse-gas supply and removal. Balance is thought to be maintained on 
habitable planets through temperature-dependent weathering reactions. Climate stability can be 
undermined by several previously-studied climate instabilities. These include atmospheric collapses 



jjlHaberle et al. 


1994; 


Read Sz Lewis 


2004) 


al. 


1997 


), greenhouse runaways 


Kasting 



albedo feedback (Roe & Baker 2010), and ocean thermohaline circulation bistability (Stommel 



1961HEPICA Community Members||2006[ ). Climate stability can also be undermined if the sign of 



the dependence of mean surface weathering rate on mean surface temperature is reversed. This 
paper identifies two new climate instabilities that involve such a reversal, and are particularly 
relevant for planets orbiting M-dwarves. 

Competition between radiative and advective heat transfer timescales sets surface tempera- 
ture on synchronously-rotating planets with an atmosphere. The atmosphere moves the surface 
temperature towards the planetary average, through radiative and turbulent heat exchange, on 
timescale t^yn- The dayside insolation gradient acts to reestablish gradients in surface tempera- 
ture, on timescale trad- We refer to planets with tdyn < Uad as dynamically equilibrated, because 
surface temperature is set by atmospheric dynamics. Venus and Titan are nearby examples. We 
refer to planets where t^yn > trad as radiatively equilibrated. Mars is a nearby example. 

Steeper horizontal temperature gradients promote atmospheric depletion if they stabilize sur- 
face liquid films, ponds or oceans in which the atmosphere can dissolve. Once dissolved, the 
atmospheric gases may be sequestered in the crust by weathering. Weathering rates are much 
faster when solvents are present and temperatures are high. Weathering and mineral formation 
can be mediated by thin films of water, and are largely irreversible on habitable-zone planets with 
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stagnant lid geodynamics (karst and oceanic dissolution layers are minor exceptions). Lithospheric 
recycling may cause metamorphic decomposition of weathering products, returning greenhouse 
gases to the atmosphere on tectonic timescales. In the absence of weathering, growth of an ocean 
can reduce atmospheric pressure through dissolution. For example, the fundamental greenhouse 
gas on Earth is CO2. The partitioning of CO2 between the atmosphere, ocean (solution) and crust 



(weathering products) is in the ratio 1:50:10^ for Earth (Sundquist & Visser 2007). Dissolution 
is fully reversible. Positive feedback occurs if reduced atmospheric pressure further steepens the 
temperature gradient. Rising maximum temperature resulting from atmosphere drawdown allow 
further expansion of liquid stability, leading to more drawdown. The zone where liquid is stable 
spreads over the substellar hemisphere. A halt to the atmospheric collapse occurs when pressure 
approaches the boiling curve, or when the liquid phase is stable over most of the dayside, or when 
thermal decomposition by crustal recycling returns weathering products to the atmosphere as fast 
as they are produced. 

Our idealized-model results motivate study of the instabilities with GCMs. 

We conclude from this study that:- 



556 1. Enhanced substellar weathering instability (ESWI) may destabilize climate on some habitable- 

557 zone planets. ESWI requires large AT^, which is most likely on planets in synchronous rota- 

558 tion. ESWI does not require strict 1:1 synchronous rotation. 

559 2. Substellar dissolution feedback (SDF) is less likely to destabilize climate. It is only possible 

560 for restrictive conditions: small oceans, highly soluble gases, and relatively thin, radiatively 

561 weak atmospheres. Furthermore, small amounts of nonsynchronous rotation can eliminate 

562 SDF. 

563 3. The proposed instabilities only work when most of the greenhouse forcing is associated with 

564 a weak greenhouse gas that also forms the majority of the atmosphere (it does not work 

565 for Earth). There are no exact solar system analogs to ESWI, although Mars comes close. 

566 Therefore, it would be incorrect to use these tentative results to argue against prioritizing 

567 M-dwarfs for transiting rocky planet searches. 

568 4. If the ESWI is widespread, we would expect to see a bimodal distribution of day-night tem- 

569 perature contrasts and thermal emission from habitable-zone rocky planets in synchronous 

570 rotation. Rocky planets with surface pressures in the unstable region would be rare, so 

571 emission temperatures would be either close to isothermal, or close to radiative equilibrium. 
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Table 1. Solar system climates, showing vulnerability to ESWI and SDF. 



Planet 


Atmospheric 
composition 


Climate-contrGlhng 
GHG 


Requirements: 
(1) 


(2) 


(3) 


(4) 


Venus 


97% CO2 


CO2 


v/ 


V 


X 


X 


Earth 


99% N2 (a) 


CO2 


X 


V 


v" 


X 


Mars 


95% CO2 


CO2 


V 






V 


Titan 


98% N2 


N2 or H2/CH4 


(b) (7) 


V 






Triton 


> 99.9% N2 


n.a. 


•y 


X 


X 


n.a. 



Note. — Requirements for the instability include; (1) The principal gats in the atmosphere 
is the main dimate-controlhng greenhouse gas. (2) The atmosphere significantly influences 
Surface temperature. (3) Removal of gas by liquid phase, via weathering or dissolution, is 
an important loss mechanism. (4) T <T^eii<-nia.x(Ts]. 

(a) Earth: Ignoring oxygen, which is an artefact of biology, (b) Titan: Small fractional or 
absolute changes in [N2] have the largest effect on Titan's r (Lorenz et al. 1999). However, 
this ignores vapor- pressure and photochemical feedbacks which may put [H2] in control 
(McKay et al. 1991). (c) Triton: Surface melting is impossible, because N^'s triple point 
is 10^ X greater than Triton's surface pressure. 
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Fig. 1. — Geometry of the idealized energy balance model for an exoplanet in 1:1 spin-orbit 
resonance. Uneven distribution of starlight (Lj,) on the planet leads to a hot (white shading, 
high Tg) dayside surface and a cool (black shading, low Tg) nightside surface. The atmosphere 
(uniform gray shading), with horizontally uniform boundary-layer temperature Tq, tends to reduce 
this temperature gradient (ATg). When Tg > Tmeit, a melt pond can form around the substellar 
point 1^ = 0, with angular radius ipmax and depth Dpond- Because rotation is slow, meridional 
winds are as fast as zonal winds, so Tg depends only on the angular distance from the substellar 
point (^). 
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Fig. 2. — Surface temperature as a function of distance from the substellar point in our energy 
balance model. Diamonds correspond to atmospheric temperature (horizontally uniform). In order 
of increasing temperature, the pressures corresponding to the diamonds are 10~^, 10^^, 10^^, 1 
and 10 bars. Arrow is the direction of increasing pressure. Radiative efficiency A = 0.1, stellar flux 
= 900 W/m^. 
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Fig. 3. — Bifurcation diagram to show the enhanced substellar weathering instabihty for radiative 
efficiency A = 0.1, stellar flux = 900 W/m^. The thick black line shows the planet-integrated 
weathering, Wt, corresponding to the temperature maps shown in Figure[2j IfVn = 0, Wt oc — (^) • 
^ = at equilibrium, and Wt equals net supply by other processes, Vn- For both the lowest and 
highest P, Wt\ as R\. Equilibria on these branches are stable. For intermediate pressures, Wt\. as 
Pt- The thick dashed line is this unstable branch. The rapid climate transitions which bound the 
hysteresis loop are shown by vertical arrows. The corresponding unstable equilibria are shown by 
open circles, and stable equilibria are shown by closed circles. The thin black lines correspond to 
Mars, Earth and Venus insolation (in order of increasing normalized weathering rate). The shape 
of the curve is explained in the text. These curves are an 8th-order polynomial fit to the model 
output. Note that as L^f) both inflection points move to higher P. 



-27- 



increased 
volcanism 



u increased 
eccentricity 

(U 
M 

a 



inward tidal 
migration 



stellar 
evolution 



atmospheric erosion, 
reduced volcanism 



substellar insolation 



50 



10- 



Radiative efficiency A=0.01 



u 
a 

0) 

u 

M 0.1 

u 



0.01 




0.001 




- Marginal climate stability 
Upper bound of liysteresis 
Lower bound of liysteresis 

- Surface water boils 

- Mean temp, exceeds 31 3K 
C02 condenses at surface 



200 400 600 800 1200 1600 2000 2600 3400 
substellar insolation (W/m^) 



Radiative efficiency A=0.3 



Radiative efficiency A=l 



50 



10 



1 - 



u 
X! 

=1 



0.01 



0.001 



m 1 
S I 

o ■ 45 




S STABLE 

8 j°---t6ES\Nr 



UNSTABLE 
TOES\N\ 




. / 



STABLE 
TO ESWI 



200 400 600 800 1200 1600 2000 2600 3400 
substellar insolation (W/m^) 



50 



10 



11 



o o. 



1 - 



a 

u 

M 0.1 



0.01 




0.001 



" I STABLE 




200 400 600 800 1200 1600 2000 2600 3400 
substellar insolation (W/m^) 



-28- 



Fig. 4. — (a) Mechanisms that can cause secular change in the location of the climate equilibrium 
Wt = Vn- (b,c,d) Habitable zone (HZ) stability diagrams for (b) A = 0.01, (c) A = 0.3, (d) 
A = 1.0. (Climates with A >> 1 are always stable against the enhanced substellar weathering 
instability). The climate states at intermediate pressure within the thick black line labelled with 
zeros are unstable to ESWI < O). Climates that approach the unstable zone from below will 
jump up to the dashed gray line. Climates that approach the unstable zone from above will jump 
down to the solid gray line. These jumps can be extreme; for example, in (b) the solid gray line 
is everywhere < 0.001 bars (and so is not visible). See the text for discussion of the speed of 
jumps. The hysteresis loop does not exist for high and and high A, and so the thick gray lines 
vanish towards the right of (d). The thin lines correspond to previously-described challenges to 
habitable-zone climate stability: moist runaway greenhouse (thin solid line) ; nightside atmospheric 
condensation of CO2 (dash-dotted line); boiling of surface water (thin dashed line). 
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Fig. 5. — Stability phase diagram, showing the effects of the enhanced substellar weathering 
instability as a function of L* and A. The jump in pressure due to ESWI can cause a runaway 
wet greenhouse (for a jump upwards in P), or a decrease in P to the triple point of water (for 
a jump downwards in P). To account for microclimates and solid-state greenhouse effects, we 
conservatively define "P < triple point" as "P < 1 mbar", which is below the boiling curves in 
Figure |4j Some curves have been smoothed with a 5th-order polynomial in order to remove small 
wiggles due to numerical artifacts. The arrow shows the change in stellar flux at lAU for a solar- 
mass star over 8 Gyr of stellar evolution in the model of Bahcall et al. (2001 ), and the circle marks 
the current solar flux. 
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Fig. 6. — Substellar dissolution feedback, for C02/seawater equilibria, A = 0.03 (note that A = 
0.03 is unrealistically low for all-C02 atmospheres), and a lOOkm-deep ocean. The vertical axis is 
P. Colored solid lines correspond to loglO(P + Ppond), i-e- the sum of the atmospheric and ocean 
inventory. Where these are equal to P, there is no ocean. Fractional ocean coverage is shown by 
the red dashed contours (contour interval is 0.1 in units of planet fractional surface area). Because 
nightside temperature is constant, fractional ocean coverage jumps from 0.5 (hemispheric ocean) to 
1.0 (global ocean). The outermost black line encloses the area where SDF is a positive feedback on 
small changes in P. Outside this area, ^^p"'' > (zero or negative climate feedback). The inner 
two contours correspond to ^^p"'' < -0.5 (strong positive feedback) and ^^p"'' < -1 (runaway). 
Runaways can only occur for deep oceans and small pond area. 



